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We show that the zero-point motion of a vortex in superconducting doped topological insulators 
leads to significant changes in the electronic spectrum at the topological phase transition in this 
system. This topological phase transition is tuned by the doping level and the corresponding effects 
are manifest in the density of states at energies which are of the order of the vortex fluctuation 
frequency. While the electronic energy gap in the spectrum generated by a stationary vortex is 
but a small fraction of the bulk superconducting gap, the vortex fluctuation frequency may be 
much larger. As a result, this quantum zero-point motion can induce a discontinuous change in the 
spectral features of the system at the topological vortex phase transition to energies which are well 
within the resolution of scanning tunneling microscopy. This discontinuous change is exclusive to 
superconducting systems in which we have a topological phase transition. Moreover, the phenomena 
studied in this work present novel effects of Magnus forces on the vortex spectrum which are not 
present in the ordinary s-wave superconductors. Finally, we demonstrate explicitly that the vortex 
in this system is equivalent to a Kitaev chain. This allows for the mapping of the vortex fluctuating 
scenario in three dimensions into similar one dimensional situations in which one may search for 
other novel signatures of topological phase transitions. 


I. INTRODUCTION 

Topologically distinct phases which cannot be classi¬ 
fied by the classical Landau paradigm comprise some of 
the most recently discovered states of matter 1-3 . An im¬ 
portant signature of these topological phases is the ap¬ 
pearance of novel, low-energy, robust, edge states; one 
such state is the so-called Majorana bound state at the 
edges of topological superconductors 4 . As ubiquitous 
signatures, the detection of these neutral fermions has 
been the main trend in the characterization of particle- 
hole symmetric topological phases. Although evidences 
of Majorana fermion physics have been identified in 
tunneling 5 and scanning tunneling microscopy (STM) 
measurements 6 , the interpretation of their signatures is 
controversial in many cases, as the imprints from the 
topological regime are often mixed with signals from dis¬ 
order and extra undesired quasiparticles. 

While the aforementioned gapless edge states act as a 
signature of topologically non-trivial regimes, the signa¬ 
tures of the transition from a topologically trivial to a 
topological phase present themselves in the bulk by the 
closing and re-opening of the excitation energy gap' -9 . 
In many of the proposed systems which can be tuned 
through a topological phase transition (TPT), the ex¬ 
citation gap is very small compared with experimental 
resolutions and cannot be probed directly. 

In this work, we show that quantum fluctuations can 
shift the spectral weight in the density of states of a 
given system before and after a TPT to further sepa¬ 
rated energies and, as a result, magnify the change of 
the spectrum resulting from this process. This situation 
will be relevant as long as the sample’s temperature is 
below HuJo/kB, where wo is the pinning frequency and 
ks is Boltzmann’s constant. We discuss this effect in 
the context of the chemical potential induced topolog¬ 


ical phase transition in the vortices of superconducting 
doped topological insulators. In this particular situation, 
we also demonstrate how the effects of Magnus forces on 
the vortex dynamics 10 have a novel signature in the spec¬ 
tral change at this TPT, exposing the pumping of vortex 
modes responsible for the phase transition, as described 
below. Our results are general, however, and can be ex¬ 
tended to other types of topological phase transitions. To 
demonstrate this, we present a way to map the 3D situ¬ 
ation into a ID setting in terms of wire networks which 
may be used to probe for the topological phase transition 
of actual Kitaev chains, Su-Schrieffer-Heeger chains and 
other unidimensional topological chains. 

To understand how quantum fluctuations affect vor¬ 
tices in superconducting doped TIs we start by discussing 
vortex dynamics in regular superconductors (SCs). This 
physics has been widely studied 11,12 and, given the nat¬ 
ural length scale of vortices, their different properties 
might display both classical and quantum phenomena. 
Within the BCS theory of superconductivity, an station¬ 
ary vortex affects the spectrum of the superconductor by 
generating in-gap modes localized around and along the 
vortex core 13 . The energy of these discrete bound states, 
known as Caroli-de Gennes-Matricon (CdG) modes, is 
given by e; = (l + |) where A is the size of the bulk 
SC gap, fi is the fermi energy and l is an integer. The sig¬ 
natures of these in-gap states have been experimentally 
observed by STM measurements 14,15 . In practice, how¬ 
ever, even though the spatial resolution of STM is well 
within the size of the vortex modes 16 , given the small size 
of their so-called mini-gap, S = the energy of each 
single mode is hard to be resolved and usually multiple 
modes are observed together 14 . 

It is well known that the pinning of vortices is neces¬ 
sary for the stability of type-II SCs. The discussion above 
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would be the final status of the problem for pinned vor¬ 
tices, were they absolutely static. Although a pinned 
vorticex has a fixed position at the sample, even at 
the lowest temperatures, their quantum zero point mo¬ 
tion cannot be ignored. Interestingly, it was shown that 
such quantum fluctuations affect the quasiparticle spec¬ 
trum, moving part of the spectral weights of the in-gap 
vortex modes to the frequencies associated with vortex 
fluctuations 17-19 . We then contend that exploiting this 
ubiquitous quantum mechanical phenomenon to probe 
for TPTs is a promising idea, leading to novel signatures 
of these transitions. 

To test this approach, superconducting doped TIs arise 
as the most natural test ground. The discovery of su¬ 
perconductivity in doped TIs triggered several studies, 
particularly because of the suggestions that doped TIs 
might realize topological superconductivity 20-24 . Theo¬ 
retical studies of superconductivity in the surface states 
of TIs started even before the experimental realization 
of bulk superconductivity in doped TIs, when it was 
shown that, theoretically, if superconductivity is induced 
in their helical surface states, vortex modes will include 
a zero-energy Majorana bound state 25 . In the context 
of bulk superconducting doped TIs, it was later shown 
that the Majorana mode at the ends of a vortex line per¬ 
sist up to a critical value of doping in these systems as 
well 26-28 . At this critical doping level, the two Majorana 
modes at the ends of the vortex hybridize and become 
gapped. The presence or absence of Majorana modes at 
the end of the vortex line contrast the two topologically 
distinct phases. In fact, the vortex in doped supercon¬ 
ducting TIs becomes effectively equivalent to a Kitaev 
chain, one of the pioneering theoretical models to realize 
topological phases and phase transitions with Majorana 
edge states 29 (check also Section V of the present work). 

As desired, the signature of this TPT also shows up 
in the spectrum of the states extended along the vor¬ 
tex. The original mechanism lies in the CdG modes. 
The important property of these states is that they are 
gapped by the small energy scale of the mentioned mini¬ 
gap. This energy protects the surface Majorana zero 
modes, confining them to the surface of the sample. 
Because of strong spin-orbit coupling and the resulting 
band inversion of TIs 30 , the Fermi surface here has non¬ 
trivial topological properties which show up as a non¬ 
zero Berry connection. The CdG modes then inherit 
this Berry phase as a modification to their energy spec¬ 
trum, which also separates in two sets due to the exis¬ 
tence of two degenerate TI Fermi surfaces, which becomes 

Ef- = ^ (l =F \ ± ) ■ Here is the Berry phase 

around the curve on the Fermi surface defined by set¬ 
ting the wave-vector along the vortex line equal to zero. 
In this case, when <!>;, = n, = 0 and the zero en¬ 
ergy surface Majorana modes at the ends of the vortex 
can merge through the gapless l = 0 mode which is now 
extended along the vortex. The richness introduced by 
spin-orbit coupling and topology in this system leads to 
the signatures that we demonstrate. 


For the physical picture of a fluctuating vortex to be 
reasonable, its position and cross-section structure must 
be well defined. Testing with some real numbers, Copper 
doped Bi^Se-i was the first topological insulator found to 
become superconducting upon doping, at 3.8 K 20 . One 
must spatially resolve the local density of states (LDOS) 
at the vortex core which, as we demonstrate, comes from 
the l = 0 and l = 1 CdG modes. Their maxima lie at 
r = 0 and are separated from the next closest mode (with 
l = — 1) by the Fermi wavevector scale r = l/fc^ ~ lOA, 
which is well within the resolution of STM. Regarding 
the energy scales, the change of spectrum at the TPT 
happens at the mini-gap energy scale S « 5 x 10 -3 K. 
This is very small compared to the spectral resolution of 
STM which is of the order of the measurement tempera¬ 
ture (3 ksT) 16 . It is then clear that an important obstacle 
to the verification of topological phase transitions in this 
system by STM is the small excitation gap. Overcoming 
this energy scale problem is the main role of the vortex 
position fluctuation we analyze. 

The paper is organized as follows. In Section II we ex¬ 
plain the model in which we base our calculations. In Sec¬ 
tion III we follow Ref. 17 , showing how the vortex fluctu¬ 
ations induce a self-energy correction which redistributes 
the peak weights in the LDOS for our specific model. 
This affects directly the tunneling conductance measured 
in STM experiments and in Section IV we demonstrate 
what are the novel consequences of this phenomenon for 
the vortex TPT in doped TIs. We believe that the ap¬ 
proach we describe in the bulk of our paper is general- 
izable to other situations and we dedicate section V to 
stipulate how to translate the ideas from the 3D context 
to ID situations concerning Kitaev chains or other linear 
or quasi-linear topological phases. We conclude in Sec¬ 
tion VI. As computations are a bit involved, we avoid 
displaying them throughout our narrative as much as we 
can. We refer the reader to the appendices, where details 
are displayed thoroughly, whenever necessary. 


II. FLUCTUATING VORTEX MODEL 


Superconductivity and the vortex quantum phase tran¬ 
sition (VQPT) in doped topological insulators may be 
understood in the weak pairing limit (£fcp 1, where £ 
is the SC coherence length) 26 . In this regime, a gradi¬ 
ent expansion can be deployed to study the effects of the 
fluctuating vortex position in the low-energy spectrum 17 . 

We start with an action of the form S = SBdG+S^ tex . 
The first term is a Bogoliubov-de Gennes (BdG) action 
for the superconducting doped TI, 

S B dG = \j d 2 rdr¥ (d T + H BdG ) * (1) 

where 


H B dG 


H T i ~ H A (r — R (t)) 
At (r - R (r)) -H T i + H 


( 2 ) 
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Figure 1. CdG vortex modes spectrum of the model from 
ref. 26 . The parameters used in the calculation are described 
in Figure 4 in the appendix. We show the energies £) + (/x) 
(blue) and E (/x) (red) for l = 0,±1 at chemical potentials 
close to /x c as 1.455. The inset displays the details for l = 0. 


To capture the coupling between electronic excitations 
and vortex fluctuations, we expand the superconducting 
pairing around the vortex rest position A (r — R (r)) ~ 
A (r) — d r A (r) • R (r). This approximation is valid at 
weak-coupling 17 , which is also the regime of validity of 
Hamiltonian (2). Within this formalism, the full problem 
is described by a perturbative action S = <S 0 + S^ tex + 
Smt- Sq is given by (1) with the BdG Hamiltonian in 
the stationary vortex limit, R (r) = 0 (explicitly given in 
(Al)). The interaction term is given by 

S^t = - j d 2 rdrR (r) • & ( ^ ) T. (5) 

The interaction between vortex modes and the fluctua¬ 
tions in the vortex position leads to a self-energy correc¬ 
tion to the energy of the CdG modes. 


Here /x is the chemical potential and the effective low 
energy 3D TI Hamiltonian is given by 

H T i = -iv D T x s.V +t z (m + eV 2 ) , (3) 

with Nambu-spinor = (ip, is y ip^ T and where ip = 
{ipA f, ipAi, ipB\, ipB\)- A, B are orbital indices and 
T-i and Si Pauli matrices act on orbital and spin 
Hilbert spaces, respectively. The superconducting pair¬ 
ing A (r — R(t)) contains a vortex profile centered at a 
fluctuating position R (r) whose dynamics is governed 
by 17 


5. 


= ^ / —Rt ( 


duj , 


eff 


27r 


iuj ) 


— LUrUJ 


UJ C UJ 

UJ 2 + CJq 


R (iuj). 

( 4 ) 


Physically, the action (4) describes a particle of mass m v 
oscillating in an harmonic trap of frequency Wo which 
depends on the properties of the trapping potential 1 '. 
This oscillator frequency dictates the qualitative features 
of the energy peak distribution of the LDOS. Finally, 
w c corresponds to a Magnus force acting on the vortex. 
The frequency u> c will be shown to play an essential role, 
introducing an energy scale for the chemical potential in 
which we have distinguished signatures of the VQPT in 
the system’s LDOS. 


III. PERTURBED LDOS 


Assuming a singlet intra-orbital pairing for doped TIs, 
the VQPT was found originally by an exact diagonaliza- 
tion of lattice toy models and a semi-classical study of 
the BdG mean-field Hamiltonian 26 , as well as numer¬ 
ically solving the self-consistent BdG equations 28 . In 
order to study the effects of vortex fluctuations on the 
LDOS, it is convenient to use a basis which diagonalizes 
the Hamiltonian at the limit of a static vortex. Thus, 
we first present the VQPT by a novel real-space diago- 
nalization of the BdG equation of Hamiltonian (Al) fol¬ 
lowing the ideas from 31 . The details follow in Appendix 
A.We expand the Grassmann fields in terms of eigen¬ 
vectors of the static-vortex BdG Hamiltonian H^ dG as 

^ = Eg=i E in Xin ( r ) ^In ( r )‘ The ei S ht arising bands 
obey H^ dG \ q ln (r) = Ef n x q lrl (r) where l and n labels con¬ 
served quantum numbers. Precisely, l represents a gener¬ 
alized angular momentum L z = —idg— , which com¬ 
mutes with the Hamiltonian (see 26 or Appendix A), while 
n labels the different eigenstates of the radial BdG equa¬ 
tion at fixed l. At weak coupling, we further project into 
the two bands which cross the doubly degenerate Fermi 
surface of Hti- Labeling these states by a = ±, we have, 
at low energies, T m xt n ( r ) ( T ) + XTn ( r ) ( T ) 

with 
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xtir, = 7sl dk 


e l ( l ^ 0 kJi- 1 (kr) 

0 

0 

e -iW ^ mk _ ^J m 2 + (/. r ) 

e“* w ^mfc + + fc 2 ) J L (kr) 


u ~L ( r ) 

V *n M 


fee l ( I+1 ) e J; + i (fcr) 


x ^ {t)= wJ ik 


C? k e lW (m k + sj m| + fc 2 ) Ji (kr) 

VK kJi -1 (fcr) 

\ 0 

/ 0 

d n e-W+WkJH .i (fcr) 

e -de (m fc - \/m| + fc 2 J J; (fcr) 

0 


V in ( r ) 


The numerical diagonalization may be done replacing 
the infinite system with a disk of finite radius R with 
a profile Ao (r) = Ao tanh (r/£) for the vortex and solv¬ 
ing the secular equation for the Fourier-Bessel coefficients 
c Ifc> d 7ki cTfc and ^"fc (details follow in Appendix A and ref¬ 
erences therein.) 

To study the VQPT we consider the lowest energy vor¬ 
tex modes. These are the CdG modes and allow fixing 
the label n —> ncd,G> which we drop. The two sectors 
(labeled by cr = ±) are connected by particle-hole (PH) 
conjugation C = p y s y IC operator (1C is the complex con¬ 
jugation operator) as C\i = XZi- The energies of the 
CdG vortex modes in this case are the expected 26 




T& (k) 

2i r 


so that E+ = —EZf Here &b(p) is the Berry phase 
calculated around the Fermi surface on the curve with 
zero wavevector along the vortex 26 . As the chemical po¬ 
tential increases, the Fermi surface enlarges and (p) 
varies from 0 to 2-7T, defining a critical chemical potential 
such that 4>b (pc) = ?r. Our results for the energies of 
the CdG modes, which are presented in Fig. 1, are con¬ 
sistent with the previous study of the phase transition in 
Refs. 26 and 28 . 

In terms of the CdG eigenstates, equation (5) is written 

S ln t = ~Y2 f dT ^i ( T ) ( T ) R ' M i/ ’ 


where 


M M' = / d2r Xi (0 f ( g r ° At ) Xf (r) • (10) 


Vortex fluctuations then generate the following self¬ 
energy for CdG vortex modes which we calculate using 
the GW approximation 32 (details follow in Appendix B), 


£f(to)= 


l' ,Q! = ± 


(iA - (agn 


Here = 


are reduced matrix elements with 


\ (M x + aiM y )^/ and = E°, + aw c /2. For 
unit vorticity, angular momentum conservation implies 
that l is connected only to l = l + al by su ch interac- 
tions. The energy scale introduced by u> v = \JZJq + w 2 /4 
( and dominated by u>q as aforementioned), represents a 
“magneto-plasma” frequency in an Einstein model 1 '. In 
Appendix B, we present closed formulas for these matrix 
elements. 

One finally needs to evaluate the LDOS, 

P(r.w) = \( e ™ M N o) 6(u)-e m ), (12) 

m,a,l 

where \N 0 ) is a A 0 -particle ground-state, \e m ) is an 
(Ao + 1 (-particle excited state (with generic quantum 
numbers rn) and ijy ; (r) is an electronic state creation 
operator at level l in sector a. Using the vortex-modes 
eigenbasis, this can be written, taking into account the 
effects of the vortex fluctuations in the self-energy, as 

fc ( r -w) = ^ (r, w) (13) 

cr=± 

fcM=--hx;— l u / (r (!!-. . (i4) 

7T U) — Ef — 2jf + It 

= £K(r)| 2 <j( W -£f-Ef(u;)). (15) 
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Through the perturbative interaction, the energy den¬ 
sity profile of CdG modes is modified with part of the 
spectral weight from u> = E° being transfered to new 
“satellite” peaks in the LDOS 17 . Both the spectrum E° 
and the profile of (r) dramatically change the phe¬ 
nomenology described by (15) when the parent metallic 
state of the superconductor comes from doped TIs, as 
compared with ordinary metals. 


IV. TUNNELING CONDUCTANCE ANALYSIS 

The local tunneling conductance is found, at low tem¬ 
peratures, by convolving the LDOS (15) with the deriva¬ 
tive of the Fermi distribution function as 

G(r,w) = ——J dto p ^r,w + u> ^ • (16) 

The normalization constant assumes an STM tip with 
constant DOS po = m e /2 tt (for a free 2D electron gas) 
with the corresponding tunneling conductance Go, and 
/ (w) is the Fermi-Dirac distribution. At very low tem¬ 
peratures, the tunneling conductance is equal to the 
LDOS, however still smoothed by the finite temperature 
effects. 

Given the atomic level resolution of STM, we can safely 
focus at the density of states at the vortex core r = 0. As 
seen in (6) and (7), the wavefunction components may 
be expanded in terms of Bessel functions. In particu¬ 
lar, at r = 0, only Bessel functions of order zero have 
non-zero amplitude while all the other Bessel functions 
vanish. From our Fourier-Bessel expansion of the CdG 
modes above, only l = 0 and, as a result of spin-orbit 
coupling, l = 1 modes have finite contributions in uj 7 (r) 
at the origin. 

The l = 0 states have energies Eq = 

cr-^ 1 + ■ These energy levels may be pumped 

from negative to positive values (and vice versa) by 


M = 


A +;a 

A 0;l 


(lo — sgn u) v — E° — uj c 


and 


£? M = 


A+'i a 

^ 1-2 


(w - sgn (^ ;,T ) 


w, 


- E% - ui c 


To find the positions of the peaks, one solves 

U-Ef- £r M = 0. (19) 

The solutions are clearly sensitive to the sign of = 
£y + aoj c /2. As we do not have an estimate for the 
actual strength of the Magnus effect, to be definite, 


changing the chemical potential, evolving the Berry 
phase from 0 to 27 t. This novel feature leads to a change 
of sign in the factors of w v in the self-energy given in 
(11) when l =0, which determine the energies of the 
satellite peaks. As a result, the TPT manifests itself 
by a discontinuous change in the density of states by 
energies of order ui v to energies of order —uj v . 

Remarkably, the local spectrum at the vortex center 
breaks particle-hole symmetry. The origin of this lies 
in the spin-orbit coupling which, together with the BdG 
doubling, filtered only the states l = 0,1 at this position, 
leaving out the l = — 1 states. Naturally the full DOS 
is PH symmetric. These points will be considered again 
in Sec. V in the context of the effective theory for the 
vortex bound states after integration in the radial and 
angular directions. 

Even more importantly, one notes that the Magnus 
force term associated with the vortex motion, whose am¬ 
plitude is proportional to to c , breaks the mirror symmetry 
which is connecting the PH sectors of the CdG modes. 
As a result, the discontinuous transition of energy of the 
CdG modes from the two er sectors does not happen si¬ 
multaneously at the same value of doping for both cases. 
This is essential for the change in the LDOS to be seen in 
this context, as it provides an energy window over which 
the density of states at the energy of vortex oscillations is 
remarkably modified by the TPT. It is also important to 
note that, for other CdG modes (such as the mode with 
l = — 1 whose maximum amplitude is at \/kp ss lOA 
away from the center of the vortex), the opposite tran¬ 
sition will happen. Given the spatial resolution of STM, 
however, the different modes should be resolvable. 

To make the above claims regarding the peak-jumping 
less abstract, let us concretely analyse the relevant con¬ 
tributions to the self-energy. As discussed, from angular 
momentum conservation, M ; + ; ; CT = 5p and from 

A”’) 7 = A^' 7 we can read the corresponding result for 
a = —. These simplifications allow us to reduce the self¬ 
energy to just a couple of relevant pieces, 


-1-;- — - 

/2) (w - sgn (-_f) oj v - E^ + u> c /2) 


/2) (w - sgn (s 0 ;<T ) lo v - E$ + u c /2) 


A 2 

we take to c = rj — . It is a simple job to notice that 
sgn (Sj t "’ <T ) = sgn (S^" ,cr ) = + and sgn (E_f) = —, for 
any value of the chemical potential. The sign of Sq ,<T , 
however, does depend on p. This allows one to define a 
value pa at which sgn (Sq ,<t ) changes. The structure of 


(17) 

(18) 
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Figure 2. Tunneling conductance for /x < p±, and /x > p± in 
blue, and p- < p < p+, in red. The large central peaks cor¬ 
respond to id « EfL 0 and id « EjL x (the energies of the CdG 
modes for stationary vortex). The inset displays the effects 
of vortex fluctuations. The smaller sattelite peaks appear at 
energies close to « ±w„. Red curves correspond to /r < p- 
or p > p + . 


Ef (ui) depends crucially on this. From the CdG spec¬ 
trum (8), we set Sq ,<t = 0 explicitly, finding 


, i I n 

-a + | a - 77 ) = 0 


<ft(/hr) = j 


■7?cr, 


( 20 ) 


where </> (p) is the Berry phase. As this phase grows 
monotonically from 0 to 27r, it is clear that the sector 
a = + has a sign change at values of p larger than those 
of the a = — sector, as long as 77 7 ^ 0. In sum, this deter¬ 
mines when each set of peaks will jump as function of p. 
In Appendix C we explore this further, also showing an¬ 
alytically that, at r = 0 , only the leftmost satellite peak 
from l = 1 will jump due to this sign change. 

Figure 2 displays our main results. It shows the dif¬ 
ferential conductivity at the vortex center G (r = 0, id) 
(more details on numerical parameters used here are 
given in Figure 4 in the appendix). Angular momen¬ 
tum conservation implies that each non-interacting en¬ 
ergy level unfolds into a set of three peaks. 

We present the differential conductance for three 
ranges of chemical potential /i < p±, in blue, and 
ft- < /x < /x+, in red, and p > p± in blue again, 
which appears to be identical to /x < p±. This hap¬ 
pens because the separation of the central peaks from 
l = 1 is E+ — E+ = S(— 1/2 + <fi(p)/2n), which can¬ 
not be resolved close to the phase transition (just as the 
peaks from l = 0 cannot be resolved at this situation.) 
In this situation, having a finite ui c is crucial to observe 
all peaks and the discontinuous effects of the topological 
phase transition. The pattern in the LDOS should be, 
for each l and sector a, of a large central peak located at 
E° with the two partners offset approximately by ifl ) 7 

with flf = y J (u> v + S + id c /2 ) 2 + 2Aj +,<T . In our case, a 
total of 12 peaks is expected for each value of the chem¬ 
ical potential (3 from 1 = 0 , another 3 from l = 1 and 


twice this due to the two sectors), not all of them being 
resolvable due to thermal effects. The large peaks closest 
to id = 0 correspond to id « E° =0 1 . The strength of the 
respective satellite peaks is suppressed by a £ -5 factor, 
where £ is the coherence length 17 . An inset displays the 
position of these peaks. 

A remarkable behavior develops in the l = 1 sattelite 
peaks (the rightmost small peaks at negative and positive 
frequencies). This is evidenced by the solitary blue peak 
at positive to. It corresponds to the contribution coming 
from id = Ep — f 2 j~, whose position jumps from this value 
by approximately 2 u> v as the chemical potential pumps 
the negative energy state at Eq into positive energies 
after crossing /x_. Similarly, when /x moves above /x + , the 
peak from id = E+ + flj 1 jumps by —2ui v . In appendix 
C we demonstrate that the approximate positions of the 
l = 1 peaks can be determined analytically. 

Concerning the magnitude of the Magnus effect, if 77 < 
1, the effects from the Magnus force are sub-dominant 
to the CdG energy gap and the sensibility to which one 
needs to tune the (zero-temperature) chemical potential 
may again be beyond technical realization at the current 
time. If 77 > 1, on the other hand, as the evolution of the 
Berry phase is from 0 to 27 t, the critical chemical poten¬ 
tials p a may not be captured as one tune /x and one will 
be bound to the regime of p- < fi < p+, which is similar 
to the standard s-wave case (except for the multiplicities 
of peaks and apparent breaking of the PH constraint.) As 
this seems to critically constrain the actual visualization 
of these effects in practice, we proceed now to consider 
some different situations in which one may actually con¬ 
trol the energy difference between the l = 1 states for 
different er = ± sectors. In this case, we will see that if 
this energy difference can be made larger, even at 77 = 0 
one may be able to capture the closing and reopening of 
the energy gap from l = 0 . 


V. ID WIRE MAPPING 


To conclude our considerations, we would like to spec¬ 
ulate about the realization of similar signatures of TPTs 
by quantum motion in other systems. Here we demon¬ 
strate concretely the claim from 26 stating that the vor¬ 
tex in superconducting doped TIs presents a topological 
phase transition equivalent to a Kitaev wire. We then 
proceed to showing that, more generally, the Hamilto¬ 
nian projected at the vortex states corresponds to a set 
of wires (or a single multiband wire) inheriting a first 
neighbor mutual coupling from the vortex fluctuations in 
3D. We then identify the important ingredients necessary 
to realize the discussed phenomena in the context of ID 
topological systems. 
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1. Vortex Hamiltonian projection 

Start with Hamiltonian (Al) from the appendix, keep¬ 
ing the z-direction terms. We also keep the vortex fluc¬ 
tuations to first order in the gradient of the supercon¬ 
ducting pairing. We have 

H BdG = H° BdG + T Z P Z - T 0 eP z 2 (21) 

-R(r)-[A-5 r A(r)] 

= H° BdG + H z +V( r). (22) 


Importantly, the sign of these couplings is the same and 
the angular integration enforces 1 = 1. For the off- 
diagonal terms we develop the couplings 

A i = f d 2 r [u+ (r)] * a z uf (r) - [v+ (r)]* a z wf (r) 

(31) 

2. ID Wire network 


We are going to project this into the lowest energy 
sectors xti ( r ) f rom (6) and (7). At finite z, we have 
XtoO) xtn cdG choosing the CdG states 

with n = n G dG- We will project the radial part the 
Hamiltonian to find out what Hamiltonian gives the 
equations of motion for ff. Considering the ± sectors 
then we have 


H ll' = Pr °J [ H BdG]u' 


(Ef 0 
l 0 Ef 


V 


H :J 


V, 

ll' 


++ 


H tu' 

H zl? 

0 

V a~ 


(23) 

(24) 

(25) 

(26) 


Notice that as (xf\p x \Xi ) = (xf\p v \Xi ) = 0, the 
fluctuating vortex potential becomes diagonal with re¬ 
spect to the ± sectors. This result is the same as we 
had found in our considerations at vanishing k z and we 
already know what this term looks like, 


K++ 

a 

0 


0 


V 


__ = R(r) • M; 


Ll' 


(27) 


with 


M+, = d 2 r 


u+ (r)^ d r Av^t (r) + v+ (r)^ d r (r)j 

(28) 


Adding up the matrix elements above gives the pro¬ 
jected Hamiltonian 


V 

R(r)- M +, 0 \ 

0 R(r)-M-J- 


( E f-efd 2 -i^id z \ 
V -iA t d z Ef + ef d 2 ) 


(32) 

(33) 


For Hermiticity M ff, = ; . For the diagonal terms 

we may still use Ef t = —Ef to write 


(Ef-e+d 2 

V -iAid z 

( E+-e+d 2 
\ -iAid z 


-iAid z \ 

E r + e i ) 

(34) 

-iAid z \ 

-(E^-efd 2 ))- 

(35) 


As the signs of ef are the same, one can easily see that 
the l = 0 Hamiltonian is essentially the same as a Kitaev 
chain. For l 0, on the other hand, the Hamiltonian 
does not describe a Kitaev chain. The PH symmetry is 
only present when both ±1, besides the a = ± sectors, 
are taken into account. In this ID projection, the con¬ 
tributions of the states in the whole radial direction are 
taken into account; in contrast, when probing the 3D sys¬ 
tem’s LDOS at the center of the vortex, we filtered the 
contributions of / = 0 and l = 1 only. Because of this, 
PH symmetry is apparently broken in Figure 2. These 
considerations are more clearly seen by writing 


M M' = 


d 2 r 


sf (r)^ d r Auf (r) + uf (r)^ d r A^vf (r) 


(29) 


Due to the vortex structure in A, we are only coupling l 
to l = / :: I. 

Now we project H z . To keep the notation short, we 
introduce 4x4 Dirac matrices a = t x u and /? = t z ctq 
as in the appendix. It is easy to see that terms linear in 
P z contribute off-diagonal in the ± sectors while terms 
quadratic with P z contribute only diagonally. For these 
diagonal terms, we develop couplings 


e -* 



ef d 2 r ju+ (r)]*/?u+ (r) - [v+ (r)]* /?v+ (rl 
ef d 2 r [uf (r)] * f3uf (r) - [vf (r )]*/3vf (r 1 )" 


Ef-jfd 2 

-iA t d z 

-iA l d z 

E f +f-T d l 

Ef — eid 2 z 

-iAid z 

-iAid z 

~ (Ef t — ef 


Ef-Efi 


PO + Pz 


e: 


Ef: 




(36) 

(37) 


(38) 


The po term does not vanish here (unless l = 0), as 
usually happens. To see that indeed the system is PH 
symmetric, one has to take into account the full second 
quantized Hamiltonian with all ±Z pairs. 



Likewise as above, the fluctuations may be written 


/R(t)-M+, 0 \ 

l 0 R (r) • M“, J 


R(r) 



Po + R (r) 


(39) 



They couple diagonally in the a = ± indices and also 
bring up an apparently PH breaking term. 

This projected Hamiltonian is then equivalent to a p- 
wave wire network. This is a very unusual network as PH 
symmetry actually connects different wires while each 
wire has PH symmetry actually broken. Although un¬ 
usual, however, similar ideas have been considered in the 
literature 34 . It is remarkable, in any case, that the 3D 
physics we started with ends up in such an exotic ID 
scenario. 


3. TPT signatures in ID 

From the above results, we may identify the minimal 
ingredients necessary to magnify the signatures of TPTs 
in ID systems by a mechanism similar to as considered 


in the vortex case. This minimal set of ingredients is 
undemanding. We list them in the context of Kitaev 
chains, as there seems to be a recent focus of interest in 
the literature concerning Kitaev wires networks (see 35 , 
for example). We stress, however, that the same ingredi¬ 
ents would suffice for other ID topological systems, such 
as Su-Schrieffer-Heeger wires or Kitaev superconducting 

1. A pair of gapped wires, one of which tunable 
through a TPT; 

2. A diagonal fluctuating coupling between them; 

With these, one may reconstruct the important features 
of Eqs.(38) and (39). This situation is illustrated in Fig¬ 
ure 3. 

Notice that the broken (PH) symmetry found in the 
projection on the vortex modes in the 3D scenario is not 
fundamental and is not included in our minimal list. It 
implies but a shift in the peaks in the spectrum, like 
as in the large peaks from l = 1 in Fig. 2, and hence 
is unimportant. Also, a single pair of Kitaev chains is 
enough (the effects from l = 2 and l = —1 in (17) and 
(18) are not important). This pair of wires could also 
be substituted by a single wire with a pair of low energy 
bands. The tight-binding model for this is written 


H l Kit ~ ~Pi ^ ° l j c 'j 
j 


l (^4M+i + a ^M4 + i + Hc ) . 


(40) 


where /q are the chemical potentials, A i are the SC pair¬ 
ings, (f>i are the corresponding SC phases and t; the hop¬ 
ping amplitude for each wire. The index l = 0, 1 labels 
the two chains. Upon BdG doubling, it is easy to demon¬ 
strate that this reduces to (38) in k-space, without the 
Po term. 

As for the fluctuation part of the Hamiltonian, one 
may have simply 

[/ = ^c° t $(r)c)+ J H.c, (41) 

j 

for a fluctuating coupling $. This should lead to similar 
self-energy corrections to the wires energies as (17) and 
(18), namely 

M = - 

LO 

and 

AT n 

S?, fc M =--, (43) 

uj-sgn [E^ J u v - E° k 

where a gives the two Nambu components. The fluctu¬ 
ation frequency uj v of $ determes the new large energy 


scale. To find the positions of the peaks, one solves again 

w - Ef ik - M = 0 (44) 

which now leads to a single sattelite peak for each energy 
level. 

Importantly, the effects of the Magnus force are not 
necessary in the ID case and, hence, a single fluctuating 
parameter is enough. This happens because one may 
(by ramping the chemical potential transversally to the 
wires, for example) keep a single wire well away from the 
phase transition with a large gap. Suppose, for example, 
wire l = 1 is kept with a large gap. In this case, the 
sattelite peaks from the two sectors in this wire will stay 
always far away from each other. This way, by tuning 
the chemical potential from wire l = 0, one can verify 
its phase transition by probing for the jumping in the 
satellite peaks of wire 1 = 1. 

As a final comment, out of the p-wave superconductiv¬ 
ity context, one might work similarly with a set of Su- 
Schrieffer-Heeger (SSH) wires. In this case, the Hamil¬ 
tonian will be similar to as the BdG Hamiltonian con¬ 
sidered so far, with the caveat that the Nambu spinor 
now should be substituted by an ordinary spinor for a 
sublattice pseudo-spin degree of freedom. The gapping 
parameters in this case will be given by staggered hop- 


A: 


0:1 


sgn 


(ur) 


- E' 


l ,fc 


(42) 
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Figure 3. ID chain minimal model to study TPTs by quantum 
fluctuations. Large circles represent lattice complex fermions 
while the prolate circles represent the Majorana fermions in¬ 
side. The l = 0 represents a Kitaev chain which we drive 
through a topological phase transition while l = 1 is kept 
at a large gap. Situations A and B correspond to the two 
deep topological regimes. The chains sites couplings are rep¬ 
resented by green wriggly lines while Majorana hoppings are 
represented by the red lines. On the right, we display an 
schematic picture of the LDOS for this system. Large black 
peaks correspond to the peaks at the energy gaps ±/q with 
k = 0, dispersion effects are not considered. As there is single 
fluctuating coupling between the chains, only a single satellite 
peak appears besides each ±/q. The purple peaks correspond 
to the satellite peaks of ±/ro. The read and green peaks differ 
by 2ui v , the fluctuation energy scale, and correspond to the 
two topological regimes A and B on the left figure. Even if 
the satellite peaks jump at the same time, since the energy 
gap from wire l = 1 is large, one may resolve the distinct situ¬ 
ations with the sattelite peaks “outside” or “inside” the peaks 
from ±/ri as in the A or B situations, respectivelly. Notice 
that PH symmetry is explicitly respected in this context. 

ping amplitudes and chemical potentials. Formally, the 
problem is the same and one may extend the results dis¬ 
cussed so far to this situation. 


VI. CONCLUSIONS 

Quantum fluctuations of vortex positions are ubiqui¬ 
tous and should manifest themselves at very low temper¬ 
atures. We found out that, in the context of doped three 
dimensional topological insulators these fluctuations may 
be exploited to magnify the signatures of topological vor¬ 
tex quantum phase transitions. This manifests at the 
LDOS at the vortex core by energy peaks which dis- 
continuously jump as function of the chemical potential. 
This finding also determined characteristic features of the 
low-energy Caroli-de Gennes-Matricon modes in this sys¬ 
tem which make them stand out as very distinct from 
standard s-wave Caroli-de Gennes modes, such as their 
spatial distribution and effects in the LDOS at the vortex 
core. Finally, our results also point to the possibility of 
capturing the effects of Magnus forces acting on the vor¬ 
tices, whose magnitude is directly related to the chemical 
potential values in which the topological phase transition 
induces peak position shifts. 

The frequency of the position fluctuations plays an im¬ 


portant role as it sets the scale of the peak jumps. In the 
context of high temperature superconductors, there are 
reports of this energy scale going up to meV 37 . It is 
important to point out that this frequency can be con¬ 
trolled to some extent and indeed increased depending 
on the properties of the vortex pinning potential. Recent 
developments in doping TIs with Niobium, which leads 
to the formation of magnetic moments in the bulk su¬ 
perconducting TI, can provide stronger pinning and so 
larger frequencies for the vortex fluctuation 33 . Measured 
physical values of the vortex fluctuation frequencies and 
Magnus force frequency in this system are not known to 
us at this point. 

Cryogenic STM measurements are fundamental to un¬ 
cover the discussed signatures. Situations with lighter 
and smaller vortices, whose zero-point motion effects 
would be stronger, could also be arranged as the vor¬ 
tex size is known to be strongly sensitive to temperature 
and magnetic field strength 36 . For vortices of too minute 
sizes, however, the Taylor expansion method deployed 
here to derive the interactions is not precise. In such 
cases, different approaches to the problem, such as used 
in 18 , are necessary in order to obtain trustworthy predic¬ 
tions. Also, a proper account for the effects of dispersion 
along the vortex may need detailed attention. It is be¬ 
yond the scope of this work to consider these. 

Finally, we studied the local physics along the vor¬ 
tex core. Projecting the Hamiltonian with the Caroli-de 
Gennes-Matricon wavefunctions we demonstrated explic¬ 
itly that the vortex line behaves as a Kitaev chain, with 
the corresponding topological phase transition. Further 
studying how the vortex position fluctuations are pro¬ 
jected into this system allowed us to find some key in¬ 
gredients which one may use to obtain new signatures 
of topological phase transitions in one-dimension. A 
promising scenario lies in the study of the density of 
states upon fluctuations of the transversal coupling be¬ 
tween a pair of neighboring gapped wires. Again, effects 
of dispersion along the wires still deserve attention. 
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Appendix A: Caroli-de Gennes modes 

Here we present our numerical method to derive 
the spectrum of vortex modes for a stationary vortex 
in doped superconducting topological insulator 26 , and 
compare the solutions with an analytical approximated 
ansatz. The latter will be used to study the effect of vor¬ 
tex quantum zero-point motion on the vortex spectrum. 
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The method which we apply is analogous to the one 
introduced by 31 and 1 ' in the context of ordinary super¬ 
conductors. We start by considering a cylinder space, 
infinite in the z direction. The Bogoliubov-de Gennes 
Hamiltonian with a static vortex centered at the origin 
reads 

H° Bd G = v D T-P-fMS + T 0 (m- eP 2 ) - A • A (r), (Al) 

where the ’’Dirac velocity” vd is set to one throughout 
our derivations and recovered to simplify the numerical 
calculations later. The Dirac matrices obey {r^Tj,} = 
2 Spv, {A a ,A b } = 2 Sab and [T M ,E] = {A a ,E} = 0. Notice 
H commutes with the kinetic Hamiltonian and is not a 
“mass” term. In our basis, a choice for the representation 
follows 

T = p z T x cr , A = pr 0 (Jo (A2) 

To = PzT z <J 0, S = PzToCTo (A3) 


Here v Pauli-matrices represent a spin-orbital coupled 
space. Noticing that 

ai = ^ d r + (A7) 

a} = - (d r - (A8) 

act as operators which lower and raise the level of Bessel 
functions (and a!a; gives the Bessel differential opera¬ 
tor itself), it is easy to find a proper basis to expand 
the states. If Ao = 0, we recover a pair of topologi¬ 
cal insulator Hamiltonians with spectra given by E )^ = 
±/i ± \J~E 2 -j-rnf, with nik = m — ek 2 and k a “radial 
linear momentum” quantum number. In the weak pair¬ 
ing approximation, since we are interested only in the 
lowest energy modes, we solve for the eigenstates of the 
TI Hamiltonian using the ladder operators above and 
project out the bands from E k + and E7~. Thus Fourier 
Bessel expanding the radial wavefunctions as 


with Nambu, orbital and spin spaces described by p, 
r and a Pauli matrices, respectively. Here A (r) = 
A 0 (r) (cos#, sin#) gives the pairing with a A 0 (r) = 
Aq tanh (r/£) profile, to be concrete. 

The vortex runs along the z direction and translation 
invariance allows us to consider the k z momentum; with 
the understanding that only k z = 0, ir are topologically 
relevant, we take k z = 0 since we are looking only for the 
low energy Caroli-de Gennes-Matricon modes 26 . 

The Hamiltonian commutes with the generalized an¬ 
gular momentum operator L z = — idg — , where 

= «Po T o <7 z = iS z - This allows writing the solution 
spinors as 


XI,n (r) = J=e-‘0-^) Vn (r), (A4) 

where l is an integer representing the standard angular 
momentum and n labels the many possible energies for a 
given l. At k z = 0, the Hamiltonian obeys a further sym¬ 
metry given by M = poT z a z . Noticing that { C,M } = 0 
and naturally {C,Hb<ig} = 0, we see that the eigenval¬ 
ues of A4 also label particle and hole partners. This al¬ 
lows one to separate </>p n (r) in four-spinors (f> l n (r), obey¬ 
ing corresponding Schrodinger’s equations with projected 
Hamiltonians H ±26 , 


Kn 


dk 


c lkh,k (r) ) 
d lk9i,k (r) J ’ 


(A9) 


where 

f f \ 1 ( hJi-i ( kr^ \ 

^’ k ^ v ( TOfc _ \/ m fc + fc2 ) ( fcr ) J 

, , 1 f (rn k + V m l + k2 ) J 1 ( kr ) ^ 

9l ' kU -Mt+Ukr) J 

(A10) 


and M kl are normalization constants given by A). = 

2 ^ k 2 + m 2 =F mfc\/m 2 + k 2 ^ / 0 °° rdrJi (kr) Ji (kr) .The 
Schrodinger’s equation reduces to 


( A +-t T+~ ) < = Kn M < (All) 

where 

Tl k , = (tm ± \J k 2 + m'l'j 6 (k - k) , (A12) 


H 




= E\ 


l,n 


£i- 


(A5) 


with respective signs, 


We focus in noticing that <j>i n = Cq n with E ln = 
-E+ The 4 x 4 reduced radial Hamiltonian reads 


H + = p z v y 


—id r + iv z 


1 

r 



m + e 



Pz + Fz 

2 ) 



- pp z - A 0 (r) p x 




(A13) 

and the spinor is 4>+ = ({cfj , {d^ k }f . 

In terms of our original variables, the wavefunctions 
are then written 


Xtn ( r ) = 


x ln 


(r) 

(r) 


+PzV: 


(A14) 
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where 


u 


+ 

In 


(r) = 


/ 


r n 

dk , Clk 

V 2 ^ + 


v 


+ 

In 


(r) = 


/ 



/ ( A . r ) \ 

0 

0 

^ e~ lW (rrik - \/m| + k 2 ^j Ji ( kr ) ^ 

(A15) 


/ e dl9 ^mfc + Vm^+F) J; (fcr) ^ 

0 

0 

y fce“*^ +1 ^ J/+i (Ax) J 

(A16) 


and the mirror (particle-hole) partners are built from 

Xin (r) = Cxtin (r). 

We then fix a finite radius R for the cylinder size which 
forces us to discretized k —> aij/R where cq 7 are the j- 
th Bessel zeroes at each l subspace. We fix a UV cutoff at 
some (large) Ao-th Bessel zero. Diagonalizing the result¬ 
ing Hamiltonian leads to the spectrum shown in 4. One 
sees two in-gap modes, one corresponding to outer edge 
modes, which we neglect, while the other corresponds to 
our desired vortex modes, as can be checked by plotting 
their respective probability densities. 

For the low energy states, n = ncdG (a label which 
we drop from now on), one may check that the spectrum 
follow the expected 


E 


± 

i 


A 2 

E f 





(A17) 


where <j> (//) is the chemical potential dependent Berry’s 
phase. At the critical chemical potential (j)(nc) = 7r - 
It grows monotonically from 0 to 2ir with the chemical 
potential. Noticing that the values of the momentum 
in fc-space are strongly localized at its Fermi value k F , 
as one might expect, it is easy to guess an analytical 
approximation for the wavefunctions which satisfies their 
desired asymptotic behaviors (see 26 for details). We have 


,+ =Ce-^ f ° dr ' A ( r ') ( f( 0 ' r \ 


(A18) 


where 

e -*(i-i)0fc F J ; _! (k F r) \ 

0 

0_ 

(' m k F - \Jrn\r +* f ) Ji ( k F r ) ) 

(A19) 

(m kF + yjml F + /t|) Ji ( k F r) ^ 

0 
0 

e -i(z+i )0k F J l+1 ( k F r) J 

(A20) 

and the new normalizations read 

N kp = 2 ( k F + m l F T \Jm 2 kF +k 2 F }j . (A21) 

Here, C is a normalization constant of order (k f/£) 7 • In 
the main text we compare the analytical and numerical 
results for the wavefunction at Z = 0 and /i ss [jl c , showing 
that the approximation indeed works. 


f(0,r) = 


^J^k f 

g{0,r ) = 


( 


\ 


o —ill 


( e 


il( 


2 ttJV, 


k F 


V 


Appendix B: Electronic Effective Interaction and 
Self-Energy 

In this section we compute explicitly the electronic self¬ 
interaction due to the interplay with the vortex fluctu¬ 
ations and the corresponding self-energy in the GW ap¬ 
proximation. 

We start from the vortex effective action of the main 
text in frequency space 


5 


vortex _ 

eff ~ 


m.y [ du , 

= - / M 


w 2 + UJq 
—U> C UJ 


LO c U) 

ur + Wg 


R (iuS) 

(Bl) 


and work at zero-temperature. Noticing that Rl (iui) = 
R (—iuj), we introduce a basis R± ( iuj ) = 
which diagonalizes the Lagrangian density as 



and the Green’s functions D F ( ico ) = 

^ (^(co ±iu> c /2) 2 + Ldlj. This sets the two impor¬ 
tant energy scales dictated by the vortex fluctuations as 
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Figure 4. Numerical results for the lowest CdG mode, l = 0, in the a = + sector. Fixing the parameters demand some care 
as one needs to consider a large enough region of k-space as to capture the TI band inversion while considering Fermi energies 
large to guarantee fc_F^ > 1 and at the same time close to the critical value /rc ~ VDy/tn/t. For all figures we use (disk 
size)/? = 5000, (number of modes) No = 300, m = 1, e = 120, Ao = 0.2, and use units with Dirac velocity vd = 15. Thus 
f = vd/ ttAo ~ 23, and the expected critical potential falls at (j,c ~ 1.37. (left) Energy spectrum as function of the chemical 
potential. A clear gap is seen at 0.2 with in gap modes. The two modes correspond to a vortex bound mode, with positive 
slope, and a gapless expected edge mode, as can be checked plotting the probability density in real space, (middle) Momentum 
space distribution of the positive slope in-gap mode at chemical potential close to the critical. The red and blue curves are 
associated with |ci = o,j| 2 and |d i=0iJ | 2 at discrete momenta j kj = ai.j/R , respectively, (right) Probability density in radial 
direction. The blue solid curve corresponds to (A14) while the dashed line corresponds to (A18), demonstrating that our ansatz 
is indeed a good approximation for the CdG modes wavefunctions. 


u) c , from the Magnus force, and u v = -\/uiq + w 2 /4, from 
the harmonic trap. 

As discussed in the former section, the low-energy 
modes divide into two Hilbert space sectors related by 


a z- mirror/particle-hole symmetry. Each sector is sub¬ 
ject to an effective potential arising after the integration 
of the vortex 0D field theory. From equations 4 we may 
write: 


e -^f,VT] oc J V[R] e- S “ +/dTE M' « (B4) 

Define U cr (iu>) = i' J V’f (iv + to) ( iv) and rewrite the scalar products in terms of the R± (ioj) coordi¬ 

nates and M a ’ a = \ (. M x + aiM y ), with a = ±. Then 


KF = 


( M ,T) 


= / d 2 


u/ (r) f <9 2 Av^ (r) + vf (r) f 5 2 A f u^ (r) 


(B5) 


This allows, with a careful consideration of positive and negative frequencies, integration over the vortex degrees of 
freedom, leading to the effective action of the electronic modes as 


dCj 


sift [r, r]=J2 (™) (*> - E n w (&) 


2tt 

cLCo 

27T 


(B6) 


i (u^D+'UZ + U^Dl'U*) 


where U" = \ (U£ + aiU y ). A tedious but straightforward simplification leads to the effective electronic self¬ 
interaction 


%[*■*]-li T, /§/£/£* 


Ll\ 


where 


Vf {iv + to) C (iv - to) Vi,i', n , n ’ (to) ipf, {iv) (iv) 

(Kf)' 


V", , {to = —— V 

l,l ,n,n v J yh / v 

" l v . 

OL — ± 


M 


uj + aiuj c /2) 2 + w 2 


(B7) 


(B8) 


From (A14), the matrix elements have a simple form 

f d 2 r [<, (r)* d g Av", (r) + < (r)* fcAt u ', (r) 


+;<r 


which also shows the convenient fact that 


MV = 

Ll' 


(B9) 
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M 7 ;CT *. 
i'i 


plane gives the self-energy of the main text 


Interaction (B8) shows a screened Coulomb-like re¬ 
tarded interaction. The self-energy in the GW approxi¬ 
mation comes now from a simple 1-loop calculation 


(«d) = -XV' mV (0) f G^M 

l' J ^ 

+E f - iw ) G T m -(bio) 

/ J Uj 


The first term vanishes. The second must be considered 
with care as the pole structure is sensitive to the structure 
of the energy levels. An integration over the complex 


l-,l' 


EE / f i— i(x‘ (j\ F<7 \ /9 

i' a=& Uw - I sgn ( I lo v +E°, J - ollo c /2 

(BH) 


where = 


M 7 

a;<T _ i,i 


and 






■ auj c / 2. 


To calculate the matrix elements one may make use of 
the Feynman-Hellman relations, adapted to our Hamilto¬ 
nian and in a finite cylinder. A long calculation making 
full use of Bessel function relations gives finally 


M + \ + 

id 


d +1 

1 ~A lj 

jj' 


(E+ +x -E+)K}+,-C l +, 


with 


-i+ij 






(B12) 


and 


JCf, =sgn(/ + l/2) (-1) 


A (Y, 


i+i a ji 


Cfj = sgn (l + 1/2) (-1) 


j-/ 2eQ A a /i+i 


i ? 3 


+ 






z+i 




'« 


a .ii\ 2 

R ) 


K 




j i+i 


(B13) 


(B14) 


Here, R is the cylinder finite radius, atji is the j-th 
zero of the 1-th Bessel function and = mj l =F 

jdsjm ?; + with nijd = m — eaji/R. The other 


matrix elements may be found from 


M a; r = - M a 't 

Id -l'-l 


(B15) 


peaks. We start rewriting, 

P(r,w) Pa (r,w) 


CT —± 


P. 


a (r,w) = --/mV 

7r z ' 


|uf (r)r 


^ w - E? - : 
z z 


ze 


(Cl) 


(C2) 


using the vortex-modes eigenbasis. STM measurements 
probe the tunneling conductance 


M «7 = 


(M+'A = fM+r) . (B16) Go r > / /\ / / /\ 

v '•'/ V / G(t,u) = --^ J du>p(r,u + u)f (to), (C3) 


These expressions are very similar to Bartosch’s, cor¬ 
rected for spin-orbit coupled states. 


Appendix C: Peak Analysis 

Here we describe in detail the determination the rel¬ 
ative sizes and positions of the tunneling conductance 


where / (w)is the Fermi distribution. 

At zero-temperature this reduces simply to the LDOS, 
up to a constant. At finite temperature we may write 

,2 


G(r, U )/0„ = -i X Y. 


4 Ml 


PO 


Z,cr=± i 


l - 


gsrK„.o) 


doj 


x/ 77,0 - w ). 


(C4) 


(C5) 
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where uj\ ^ q is the i-th solution to 

u-Ef- Sr (w) = 0. (C6) 

This represents a cubic equation, thus with three so¬ 
lutions. While (C6) determines where are the relative 
positions of the peaks in energy space, the derivatives 

dY ‘ l will fix the peaks relative sizes. 

We focus most of our analysis at |r| =0, which, from 
(A14), means that only the states with l = 0, 1 give non¬ 


vanishing contributions. The relevan self-energy contri¬ 
butions were considered in the main text in equations 
(17) and (18). To determine the relative sizes and po¬ 
sitions of the peaks, we examine the derivatives of the 
self-energy, as well as equation (C6) explicitly. 

1. Peak sizes 

The derivatives of the self-energies read, after some 
simplification 


(w) 


/4 +;<T 

^1,2 

A +;a 

A 0;l 

(C7) 

dto 

(Awf 

- 8 - w c /2 - w„) 2 

(Aw) 7 + 8 + w c /2 — asgn (/i — p, a ) w„) 2 

dX z M _ 


A +;a 

^0:1 

A 

+ 

-6^ 

(C8) 

doj 

(AwJ 

- 8 - w c /2 - w„) 2 

(Acjq + S + cu c /2 + lu v ) 


where A uif = to — E° and 8 is the mini-gap. 

The matrix elements are much smaller than the other 
physical quantities. Dimensional analysis and explicit 
manipulation of (B5) shows that, at constant u) v / Ao, 
these overlaps sizes depend on the coherence length as 
£ -517 . The peak sizes, nevertheless, are going to be 
sensitive to A“T. As will be seen in the next subsec¬ 
tion, the satellite peaks positions are dominated by the 
vortex oscillation frequency u> v . Plugging A uf ~ 0 or 
A uf « ±uj v one sees that cffif (oj) /dco is small (con¬ 
cretely it isoc < 1) at Aw^ « 0 while it may 

be larger at Acjg ss ±uj v , going as ~ —Aj) 4 ' [^], where 

s = S+ 2 u^ 2 ■ The latter case reduces greatly the size of 
the satellite peaks from l = 0, similarly as pointed by 
Bartosch et al. 17 . 


2. peak positions 

Our last goal is to explain the positions of the peaks 
as function of the chemical potential, demonstrating that 
they are much less sensitive to the matrix elements than 
the peak sizes and that they are mainly fixed by the vor¬ 
tex fluctuation frequency, which might be much larger 
than the other energy scales of the problem. 

Simplifying the self-energy and plugging into (C6), 
shows that independent of chemical potential, for l = 0 
we have 

Awq (Awq ) — (<5 + ui c /2 + u! v ) + (Aq.’j 4 + 

+ (<W + 5 + u) c / 2) (Ag.’-j 4 — jt 0 ) = 0. 

(C9) 

Using Aq.’j 4 « 0 we get results similar to reference 17 

for ordinary s-wave superconductor. Since the matrix 


elements are much smaller than the other parameters, 
we can neglect them in above equation. We then get 

AixIq (Awq ) — (5 + u) c /2 + u> v ) — 0, (CIO) 

for any ^. This gives a central and two satellite peaks at, 
respectively 


AWg = 0 

(Cll) 

Awq = (u! v + 8 + U! c / 2) 

(C12) 

Awq = — (c o v + 8 + u> c /2). 

(CIS) 

For 1= L, we may as well neglect the contributions 

from the matrix elements. For /i < /r_, 

Auq (Awj — co v ) 2 — (8 + w c /2) 2 

= 0 (C14) 

Aw 4 (Aw) 4 )” — (u! v + 8 + w c /2) 2 

= 0. (C15) 

So we have peaks at 

Auq = 0 

(C16) 

Auq — cu v T (^ A cu c /2) 

(C17) 

A u>i = w v — (5 + ui c /2) 

(C18) 

and 

= 0 

(C19) 

Ac — (c u v S -|- cj c /2) 

(C20) 

Ac— — iyiy H - S (^d 2) • 

(C21) 


For n > p> + the role of + and — in above equations are 
exchanged. Since the total density is the sum of contri¬ 
butions form both a = ± sectors, and the gap between 
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E+ and Ep goes as <5(1/2 — <j>/2ir), the LDOS in the two 
regimes of < fi- and /r > p+ look the same. 

We now get to the most important regime of /!_ < p < 
p + . The position of the peaks for both a = ± sectors are 
at 


A w? = 0 (C22) 

A.wf = (tu v + 5 + w c /2) (C23) 

A.wf = - (w„ + <5 + w c /2). (C24) 

Clearly, as /i crossed p- the third peak for <r = — sector 


is shifted by —2cu v and this leads to a clear modification 
of LDOS which persists up the /i = p + at which the 
peak form the a = + sector moves by 2co v and recovers 
the original LDOS. 


The “creation” of a satellite peak at positive energy 
should not happen without an accompanying compen¬ 
sation of a positive energy peak jumping into negative 
energies. Indeed, such a compensation does occur for the 
contribution of l = — 1 (which exchanging angular mo¬ 
mentum with the vortex motion is connected to l = —2 
and l = 0, the latter giving the jump.) It just turns out 
that, since the spatial dependence of the LDOS is de¬ 
termined by '«/ (r), as can be seen from (C4), the peaks 
from l = — 1 do not contribute to the LDOS at the center 
of the vortex, r = 0. The peaks from l = — 1 should con¬ 
tribute to the LDOS at a distance ~ fcjj 1 from the vortex 
center, which should be of the order of ten Angstroms 
in a superconducting TI. This can be resolved with the 
current STM technology. 
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